Temperature dependent fluorescence in disordered Frenkel chains: interplay of 
equilibration and local band-edge level structure 
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We model the optical dynamics in linear Frenkel exciton systems governed by scattering on static 
disorder and lattice vibrations, and calculate the temperature dependent fluorescence spectrum and 
lifetime. The fluorescence Stokes shift shows a nonmonotonic behavior with temperature, which 
derives from the interplay of the local band-edge level structure and thermal equilibration. The 
model yields excellent fits to experiments performed on linear dye aggregates. 
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Optical dynamics and emission of light in solid state 
systems are topics of considerable fundamental and tech- 
nological importance. A variety of materials is currently 
studied as light emitting devices. Examples are semicon- 
ductor quantum wells and wiresjlj, conjugated polymers 
0, and molecular aggregates Q. Closely related is the 
quest for understanding the optical dynamics in photo- 
synthetic antenna systems 0, y| . The fundamental un- 
derstanding and description of the emission properties of 
all these systems are challenging, because they involve a 
complicated interplay of processes. The emission is usu- 
ally caused by excitons and is particularly sensitive to 
the scattering of these quasi-particles on static disorder, 
which may lead to their localization, and lattice vibra- 
tions. The latter results in vibronic relaxation and hop- 
ping of the excitation between (localized) exciton states, 
which gives rise to spatial and spectral diffusion. These 
diffusion processes compete with spontaneous emission 
and nonradiative loss. Depending on the details of this 
competition, the excitons may or may not reach thermal 
equilibrium before emission occurs. Temperature plays 
an important role in this competition. The thermal de- 
struction of cooperative emission in molecular aggregates 
, quantum wells , and quantum wires 0, is a well- 
known example of an interesting temperature dependent 
effect. The Stokes shift of the exciton fluorescence may 
exhibit interesting thermal behavior as well 0, |nj . 

In this Letter, we theoretically study the exciton dy- 
namics and the resulting emission properties of linear 
chains carrying small-radius excitons. We account for 
both exciton localization by static disorder and scatter- 
ing of the excitons on lattice vibrations. The interplay of 
nonequilibrium conditions and the local band-edge level 
structure of disordered chains leads to interesting effects 
in the emission. A comparison to experiment is made. 

We consider a chain of N molecules described by the 
Frenkel exciton Hamiltonian 
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where e„ is the excitation energy of the nth molecule and 
|n) denotes the state in which molecule n is excited and 
all others are in their respective ground states. In order 
to model the effect of random solvent shifts, the e„ are 
assumed to be mutually uncorrelated Gaussian stochastic 
variables with mean eo and standard deviation a. The 
hopping integrals J nm are nonrandom and originate from 
dipolar coupling: J nm — —J/\n — m| 3 , J nn = 0. Here, 
J is positive, implying that the optically allowed exciton 
states reside at the bottom of the exciton band. 

We will assume that the coupling of the excitons to 
lattice vibrations is weak. It is then appropriate to use a 
basis of eigenstates of H. This basis consists of N exciton 
states, which are labeled by v = 1, . . . , N. The energy of 
the vth exciton, E v , is the vth eigenvalue of the N x N 
matrix H nm — (n\H\m), while the nth component of the 
corresponding eigenvector, tp^n, gives the amplitude of 
the vth exciton state on molecule n: \v) — Y] 1 yw|n). 
The disorder a leads to localization of the eigenstates and 
gives rise to a tail in the density of states (DOS) below 
the bare band edge 0, 0. This tail is of particular 
interest for optical properties, as the states residing in it 
carry most of the oscillator strength, F v = \^2 n fvn\ 2 - 
While the DOS accumulated over the entire chain is a 
smooth function of energy, locally it exhibits structure 
(l4|. This local level structure is of key importance to 
the low-temperature dynamics and is described below. 

The exciton states in the tail of the DOS are localized 
on finite segments, which near the band edge have a typ- 
ical size N* (see Fig. QJ. Within its own segment, each 
one of these states is the lowest-energy eigenstate. They 
have wave functions without nodes in the segments and 
carry an oscillator strength F v « N* . Their spontaneous 
emission rates, 7^ = joF v (with 70 the emission rate of 
a single molecule), thus scale like N* . This is referred 
to as superradiant enhancement 0. Bottom states on 
neighboring segments have small overlap. Often a sec- 
ond, higher- energy, exciton state may be distinguished 
on a segment, separated from the bottom state by an en- 
ergy of the order of 3tt 2 J/N* 2 , as is dictated by energy 
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FIG. 1: Schematic picture of the localized exciton wave func- 
tions near the lower band edge in a disordered linear chain. 
Also indicated are the optical pumping and emission pro- 
cesses (wavy arrows) and various intraband relaxation pro- 
cesses (dashed and dash-dotted arrows). 

quantization. These states have one node within the seg- 
ments and small oscillator strengths. Still higher local 
states are usually delocalized over two or more segments 
and carry negligible oscillator strengths. The local level 
structure is washed out in the average over the entire 
chain due to the fact that the distribution of energies of 
the bottom segment states also has a width of the order 
of 3ir 2 J/N* 2 Q. 

Weak exciton-phonon coupling leads to scattering be- 
tween the exciton eigenstates. To describe this dynamics, 
we will use the Pauli master equation for the populations, 
P v {t), of the localized exciton states: 

N 

P v = R v {t) - lv P v + £(W„ M P M - W^P V ) . (2) 

Here, the dot denotes the time derivative, i?„(i) is the op- 
tical pumping rate, j v is the above denned spontaneous 
emission rate, and W v>1 is the transition rate from the 
localized state \i to the state v induced by the exciton- 
phonon scattering. For the latter, we will use 

N 

W v , = Wo S(E U - E„) G(E V -E^Y, vlnVln > ( 3 ) 

n=l 

where the constant Wo characterizes the scattering, the 
sum over sites is the overlap integral of exciton proba- 
bilities for the states fj, and u, and S(E V — E^) describes 
that part of the E v — E^ | dependence of W v ^ that derives 
from the energy dependence of the exciton-phonon cou- 
pling constant and the phonon density of states. Finally, 
G(E) = n(E) if E > and G(E) = 1 + n(-E) if E < 0, 
with n(f2) — [exp(f2//cBT) — l]- 1 the mean occupation 
number of a phonon state with the energy Q. Equation 
©, which meets the principle of detailed balance, is the 
generic result from first-ord er p erturbation theory in the 
exciton-phonon interaction LLa, Ljj • 

We have numerically simulated the above model of 
intra-band exciton dynamics for scattering on acous- 
tic phonons of the chains' host medium. We have set 



S = \E V — Efj\/J, which properly accounts for a decrease 
of the coupling of excitons to acoustic phonons at small 
wave vectors and prevents the divergence of W v ^ at small 
energy differences. We will illustrate the salient results 
of our model by using parameters that are relevant to 
molecular aggregates of the dye THIATS: eo = 18138 
cm- 1 , J = 740 cm" 1 , and 7o = (2/31) x 10 s s" 1 [H03. 
These parameters are fixed; the two free parameters are 
a and Wo, for which we will use the values that yield the 
best fit to experimental data: a = 0.23J and Wo = 30J, 
respectively (vide infra). A more elaborate analysis of 
our model, covering a wider rang e of values for a and Wo, 
will be discussed elsewhere .19]. The motivation to fo- 
cus on the example of THIATS aggregates, is that a large 
number of temperature dependent spectroscopic data are 
known for these systems [ill l2p| . We will limit our at- 
tention to temperatures small compared to the Davydov 
splitting of 3000 cm- 1 that exists for these aggregates 
between the so-called J band at the lower exciton band 
edg e and the so-called H band at the upper band edge 
[111 llS| . For these temperatures our effective model with 
one molecule per unit cell may be used (SJ. We have 
used a chain size of N — 500 in our simulations, which 
far exceeds the typical band-edge exciton size N* = 30 
for a = 0.23J. 

We will first address the absorption and steady-state 
emission spectra, which are calculated through 

1 N 

ME) = ^(Y-y^E - E «)) ■ ( 4a ) 

respectively. Here, the angular brackets denote the av- 
erage over the disorder realizations. Furthermore, P^ is 
the steady-state solution of the master equation Eq. (J2J 
with the source term given by P v (t) = F v if E v falls in- 
side a narrow window in the blue wing of the J band and 
R v = otherwise. This blue-tail excitation condition 
agrees with common experimental conditions. 

The spectra calculated for T = 10 K are shown by 
solid lines in Fig. |3 together with the J band and fluo- 
rescence spectrum measured for THIATS aggregates at 
the same temperature (dots) ^lj- Clearly, both experi- 
mental spectra are fitted very well within our model. The 
main deviations occur in the tails, which are quite sensi- 
tive to the precise nature of the disorder (diagonal versus 
off-diagonal disorder, no perfect alignment of molecular 
dipoles). We note that the fluorescence band is shifted 
to the red (Stokes shift) and narrowed with respect to 
the J band. Keeping in mind that the absorption band is 
dominated by the set of superradiant bottom states of the 
various localization segments on the chain, this is an indi- 
cation that after the fast initial relaxation to a particular 
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FIG. 2: Calculated absorption and fluorescence spectra (solid 
lines) at 10 K for the model parameters described in the text, 
compared to the J band (o) and fluorescence spectrum (•) 
for THIATS aggregates measured at T = 10 K |Ti|. The 
dashed line is the fluorescence spectrum calculated for the 
same parameter set, except that now Wo = 10 -3 J. 

bottom state (e.g., the level a in Fig.^), the excitons have 
a chance to relax to even lower states. This is a signature 
of migration to other localization segments, a process 
which at low temperatures is hindered by the small over- 
lap between exciton wave functions of different segments. 
The fact that, in spite of this, migration may occur before 
radiative emission takes place, is due to the large value of 
Wo, which brings the system into a regime that we refer 
to as the fast-intersegment- relaxation limit. To illustrate 
the role of Wo, we also present in Fig. [2] (dashed line) 
the simulated fluorescence spectrum for Wq — 10~ 3 J (all 
other parameters as above) . In this slow-relaxation limit 
the fluorescence spectrum closely follows the absorption 
band and no Stokes shift occurs. 

At zero temperature, the possibility for the excitation 
to migrate to other localization segments, after the initial 
relaxation process, is very small. The reason is that the 
migration may then only occur via a jump to a lower ex- 
citon state on a neighboring segment, for which at least 
some (albeit small) wave function overlap exists (cf. the 
process denoted by the dashed arrow in Fig. ^| . How- 
ever, the probability of finding a lower exciton state on 
a neighboring segment is small, because the initial state 
already lies in the tail of the DOS. Therefore, in practice, 
low-temperature exciton migration is limited to at most 
one jump. As a consequence of this blockade for diffu- 
sion, there is not much opportunity for the excitation to 
explore the full DOS during its lifetime and to equilibrate 
to the lowest state on the disordered chain. Thus, even 
for large Wq, the zero temperature Stokes shift is limited 
to values of the order of the J band width. 

Upon increasing the temperature, higher states may 
thermally be excited and one expects the Stokes shift to 
decrease and eventually vanish. This is the generic be- 
havior for single molecules. Our simulations show that 
at elevated temperatures this also occurs for the exci- 
ton system (Fig. EJa)). However, our results also re- 
veal a counter-intuitive intermediate regime, where the 
Stokes shift increases with temperature. This nonmono- 
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FIG. 3: (a) Calculated Stokes shift of the fluorescence spec- 
trum as a function of temperature for the parameter set de- 
scribed in the text (•). The solid line is a guide to the eye. 
Diamonds mark the measured Stokes shift for THIATS aggre- 
gates reported in Ref. [Tl|. (b) As in (a), for the fluoresecence 
lifetime instead of the Stokes shift. Diamonds are the exper- 
imental data from Ref. [20|. 



tonic behavior is intimately related to the interplay of 
equilibration and the local level structure near the exci- 
ton band edge. Upon slightly increasing the temperature 
from zero, a small fraction of the population is thermally 
raised from the bottom state of a segment to a higher 
exciton state that overlaps with the same segment. As 
this higher state generally overlaps with more segments, 
this provides a means of distributing the population over 
two or more neighboring segments. The total process 
is denoted by the dash-dotted arrows in Fig. ^ This 
overlap-assisted activated migration is slow, but it allows 
the excitation to probe more segments for a low-energy 
position and circumvents to some extent the blockade for 
migration which occurs at T — K. Hence, the excitons 
get closer to thermal equilibrium, which at low temper- 
atures leads to a larger Stokes shift. Such anomalous 
low-temperature behavior has also been predicted and 
observed for Wannier excitons migrating between local 
potential minima in narrow quantum wells |lfj| . 

As we also see in Fig.l^a), our numerical results, ob- 
tained for the same parameter set as was used to fit the 
low-temperature spectra, agree very well with the experi- 
mental data for THIATS aggregates reported in Ref. • 
This gives strong evidence for the validity of our dynamic 
model and the concept of overlap-assisted thermally- 
activated intersegment hopping. Using the typical local 
level separation 37r 2 J/iV* 2 and the superradiant emis- 
sion rate N*jq, it may be shown that the temperature 
interval over which the Stokes shift behaves in an anoma- 
lous way (i.e., increases with T) scales roughly linearly 
with the absorption line width and only weakly (inverse 
logarithmically) depends on the ratio Wo/70 |l9|. 

We finally turn to the lifetime of the total fluorescence 
intensity after pulsed excitation. For the generally non- 
exponential intensity trace I(t) = —(^2 v Pu{t)), we de- 
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fine the lifetime as its first moment, i.e., the expectation 
value of the photon emission time: r = f P v (t))dt 
[Tt| . Here, P v (t) is the solution of Eq. J2J with R u {t) = 
cF v 8{t) if E v falls in a narrow window in the blue tail 
of the absorption band and R v {t) = otherwise. The 
normalization constant c is such that directly following 
the excitation pulse the total population on the chain, 
^2 V P V , equals unity. 

The generic behavior of the lifetime as a function of 
temperature is well known. At low T, the decay is domi- 
nated by the superradiant bottom states of the localiza- 
tion segments, leading to a fluorescence lifetime of the or- 
der of l/N*"f . With increasing T, the exciton population 
is shared by the superradiant states as well as an increas- 
ing number of higher- lying dark states. This leads to an 
increase of the fluorescence lifetime, as has been observed 
for various exciton systems H 0, ll H El E2 . For THI- 
ATS aggregates the lifetime shows a nearly linear tem- 
perature dependence from to 130 K j2^| (see Fig.^b), 
diamonds). For BIC aggregates a similar linear increase 
was observed and, based on the flat nature of the DOS 
for ideal two-dimensional systems, this was put forward 
as evidence for a two-dimensional structure of these ag- 
gregates [i^. It turns out, however, that the interplay of 
the (local) DOS and equilibration in our one-dimensional 
model with disorder may equally well explain the quasi- 
linear behavior. The solid line in Fig. ISIb) gives our 
simulated results, obtained for the same parameters that 
were used to fit the low-temperature spectra of THIATS. 
We conclude that the model, without any new parame- 
ters, yields an excellent fit to the experimental data. We 
also mention that at small Wo (slow relaxation), the lin- 
ear temperature dependence may be preceded by a low- 
temperature plateau [l9| , which is the type of behavior 
that is observed for aggregates of pseudo-isocyanine . 

In summary, we have analyzed the temperature depen- 
dence of the fluorescence in one-dimensional disordered 
exciton systems and obtained excellent fits to experimen- 
tal data for aggregates of the THIATS dye molecules. As 
far as we are aware, this is the first model calculation that 
yields such good agreement with a variety of fluorescence 
experiments on linear dye aggregates. Previous calcula- 
tions focused only on the temperature dependent lifetime 
and did not account for exciton localization by disorder 
H. Our approach, based on the master equation with 
perturbatively calculated rates W UI1 , is consistent as long 
as the rates for scattering between the optically dominant 
states do not exceed their energy separation. For the pa- 
rameters used for THIATS we find that, in spite of the 
seemingly large value of Wo, at low temperature this cri- 
terion is just obeyed, due to the small overlap integral 
Y2n=i finf'fm ~ 1/A* for the exciton wave functions 
near the band edge. With increasing temperature, the 
scattering rates increase and around 100 K they lead to 
a visible contribution to the absorption line width (ho- 
mogeneous broadening) ^lj. At such high temperatures 



the master equation description breaks down and the full 
exciton density matrix must be considered instead Q. 
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